Liesegang patterns: studies on the width law 
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The so-called width law for Liesegang patterns, which states that the positions x n and widths 
w n of bands verify the relation x n ~ iu™ for some a > 0, is investigated both experimentally and 
theoretically. We provide experimental data exhibiting good evidence for values of a close to 1. The 
value a = 1 is supported by theoretical arguments based on a generic model of reaction-diffusion. 
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I. INTRODUCTION 

In recent years, pattern-forming chemical, physical and 
biological processes came to the foreground. The three 
main causes of this strong interest are the exotic nature 
of these processes, designed pattern-forming experiments 
turned up in physical and chemical laboratories and, last 
but not least, the dramatic evolution of computing power, 
allowing to perform extensive numerical simulations as 
well as to solve numerically the reaction-diffusion differ- 
ential equations which underly these phenomena [EJ. 

When an electrolyte diffuses into a gel containing an- 
other electrolyte, the eventual formation of a rhythmic 
pattern of precipitate (a weakly soluble salt) is known as 
the Liesegang phenomenon J^||. This pattern strongly 
depends on the dimensionality as well as on the geometry 
of the experimental setup. 

The simplest and usual way to perform a Liesegang 
experiment is to take a tube and fill it with a gel which 
contains the so-called internal electrolyte dissolved in it. 
The other reactant, refered to as the outer electrolyte, 
is poured onto the surface of the gel at one extremity of 
the tube. Usually, the concentration of the latter is much 
higher than that of the inner one. The main role of the gel 
in this experiment is to prevent convection of solutions 
and sedimentation of precipitate. Thus the only mech- 
anism allowed here for the transportation of material is 
diffusion, i.e. the whole process is diffusion-limited. 

Although the Liesegang phenomenon has been known 
for more than one hundred years, there are still many dis- 
cussions on the detailed mechanisms responsible for these 
static structures. Nevertheless, most of the Liesegang 
patterns one obtains in experiments or observes in na- 
ture are characterizable by the following generic laws: 

1. The ratio between the positions x n , x n+ ± of two 
successive bands, measured from the surface of the 
gel, tends to converge to a finite value that is usu- 
ally larger than one (spacing law) B: 



2. The ratio between the position of a band and the 
square root of the time elapsed until it appears, 
denoted t n , is constant. This law is a direct conse- 
quence of the diffusive nature of the dynamics (time 
law) @: 

3*71 ^ ^/tn (2) 



3. It is also generally observed in experiments that the 
width w n of the bands increases with their position 
as a power law (width law) raj: 



(3) 



with a > 0. 



x n+l n large 



Until now, a lot of attention has been paid to the spac- 
ing law (see §,[§,§, @). On the contrary, the width 
law has largely been ignored for several reasons. Fom the 
experimental point of view, it is difficult to make precise 
measurements of the width of the bands without hav- 
ing recourse to sophisticated digitizing methods. From a 
theoretical point of view, a complete description of the 
dynamics of bands formation requires the knowledge of 
the detailed mechanisms involved in the coarsening pro- 
cess. Accordingly, two different predictions for a have 
been previously obtained fl|lT|. 

It is our purpose, in this paper, to focus on the width- 
law aspect of the Liesegang phenomenon, and to try 
hereby to remove some ambiguities that still surround 
this basic feature of this class of systems. In section II, 
we will expose and discuss the results of Liesegang ex- 
periments performed by one of us (M.Z.), for which the 
widths of the bands were measured. The analysis of the 
data provides an accurate estimation of the exponent a. 
In section III, a theoretical argument, independent of the 
nature of the detailed mechanisms involved in band for- 
mation, is developed to support the experimental find- 
ings. Finally, concluding remarks are made in section 
IV. 



l+p 



(1) 



with p > 0. 



II. EXPERIMENTAL PART 
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A. Materials and Sample preparation 

We have experimentally studied the formation of 
Mn(0H)2 precipitate by diffusion of NH4OH into chem- 
ically cross-linked polyvinyl(alcohol) (PVA) hydro-gel 



containing the MnSC>4 electrolyte. The following reac- 
tion occurs when the spatially separated reactants come 
into contact: 

Mn 2+ + 20H- # Mn(OH) 2 (4) 



FIG. 1. Picture of the resulting pattern obtained in one of the three experiments described in Section 2. 



The highly swollen PVA-hydrogels were prepared by 
cross-linking of primary PVA -chains with glutaric alde- 
hyde (GDA) in aqueous solution, containing dissolved 
MnS0 4 . Commercial PVA (Merck 821038) and solu- 
tion of 25mass% GDA (Merck) were used for prepara- 
tion. The initial polymer concentration, as well as the 
cross-linking density, were kept constant. The polymer 
content of the PVA gels was in every case 4.0 mass%. 
The ratio of monomer unit of PVA (VA) to the cross- 



linking agent, GDA was: (VA)/(GDA)= 250. In order 
to study the precipitation and band formation in swollen 
networks, one of the reactants, MnSO/t, was mixed with 
the polymer solution containing the cross-linking agent. 
The concentration of MnS04 in the gel phase was 0.03 
mol/dm 3 . The gelation was induced by decreasing the 
pH of the system by H2SO4 (Carlo Erba). The solution 
was then poured into glass tubes. The gelling process 
usually took 3-4 hours. 



Sample label 


ao [mol/dm?] 


60 [mol/dm^] 


a 


A 


3.0 


0.03 


0.94 


B 


5.0 


0.03 


0.99 


C 


14.7 


0.03 


0.90 



Table I. Concentrations used in each experiment, ao is the concentration of the the outer electrolyte (NH4OH) 
and 6q, the one if the the inner electrolyte {MnSOa). The last column gives the width exponent a. 



The tubes were sealed, placed vertically, and allowed 
to stand undisturbed at a constant temperatures of 17°C. 
After completion of the network formation, the gels were 
brought into contact with the outer reactant, NH4OH, 
thus defining the initial time t — for the experiment. 
The concentration of outer electrolyte was varied between 
3 and 14.7 mol/dm 3 . We maintained the constancy of 
boundary condition by refreshing the outer solution con- 
tinuously. Let us underline that due to the choice of the 
initial conditions, these experiments are effectively one- 
dimensional. The relevant space coordinate is the one 
that goes along the longitudinal axis of the tube. 



a real-time video digitizer card. We also used a scanner 
to digitize the experimental results. 

We first determined the coordinate of the gel surface. 
Then, the positions of the N bands, {x n } n =i,.N, were 
evaluated. Note that x n stands here for the position of 
the upper limit of the precipitate zone. We finally mea- 
sured the positions {x n }n=i..N of the lower boundaries of 
the precipitate layers. The difference of these two quanti- 
ties provides the width of the n-th band : w n = x n — x n . 

It must be mentioned that due to the sharp boundaries 
of the precipitate regions, the uncertainties on the width 
values are estimated ranging from 2% to 5%. 



B. Experimental data analysis 

The penetration of NH4OH into the gel containing 
MnS04 results in a sharp boundary interface of Mn(OH)2 
precipitate bands. In many cases we obtained more than 
20 bands. Figure 1 shows a typical example. 

In order to determine the space coordinate of the bands 
obtained, we used a digital video system. A CCD camera 
with an 1/3" video chip was connected to a PC through 



C. Experimental results 

Figures 2-4 show the results of three different Liesegang 
experiments for which the concentrations of the elec- 
trolytes are given in Table 1. They take the form of 
log-log plots relating the widths w n to the positions x n 
of the bands, as measured from the digitized images. The 
dotted lines are least square linear fits whose slope gives 
the width exponent a. These slope values are shown in 
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Table 1. For all the samples, the width exponent is near a = 1. 



A 5.5 - 











y' 

,i 

y* 


- 








*-* 








a = 0.99 




¥' 







9.5 
ln(x„) 



FIG. 2. Log-log plot of the widths w n of the bands versus their positions x„ for the sample A. Lengths are in arbitrary 
units. The dotted line corresponds to a least square linear fit. 
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FIG. 3. Log-log plot of the widths w n of the bands versus their positions x n for the sample B. Lengths are in arbitrary units. 
The dotted line corresponds to a least square linear fit. 
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FIG. 4. Log-log plot of the widths w„ of the bands versus their positions x n for the sample C. Lengths are in arbitrary units. 
The dotted line corresponds to a least square linear fit. 
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In these experiments, the two chemical species react 
and produce an intermediate compound (called C) which 
also diffuses inside the gel. Its associated traveling front 
can be put in evidence by light transmission techniques. 



III. THEORETICAL APPROACHES 

A. Two classes of competing theories for Liesegang 
patterns 

The several mechanisms one can consider for explain- 
ing the formation of Liesegang patterns will be, in the 
context of this paper, classified in two categories p0[ . 
Both have in common the following set of "free" param- 
eters : 

• i). The concentration of the outer reactant A, de- 
noted do- It is generally kept fixed at the open 
extremity of the tube during the whole duration of 
the experiment. It therefore leads to the follow- 
ing initial- and-boundary condition for the field a 
describing the concentration of species A : 

o(0 < x < L,t = 0) = (5) 
a(x = 0,t > 0) = a (6) 

where L stands here for the length of the test-tube 
hosting the gel. 

• ii).The initial concentration of the inner electrolyte, 
60. One has: 

b(0 < x < L,t = 0) = 6 (7) 

• iii) . The diffusion coefficients of both species in the 
gel, D a and D b . 

Each of these two classes of models we distinguish here 
try to capture a certain family of microscopical mecha- 
nisms leading to the formation of bands. In this respect, 
they both can be characterized by their own set of pa- 
rameters associated to these particular mechanisms : 

1. In the first category, we put the ion-product su- 
per 'saturation type theories. They have in com- 
mon the fact that the outer and inner reactants 
directly turn into precipitate where their local con- 
centration product reaches some threshold q* . One 
can symbolize this process by the reaction scheme 
A + B — > D_ where D_ stands here for the precipitate 
state. The typical and most studied scheme of this 
type is Prager's theory for periodic precipitation 
H , firstly proposed by Wagner M , and investigated 
later also by Zel'dovitch et al M. 



2. The second category contains what we will call the 
intermediate species-type theories. We group into 
this set all the schemes that involve the existence 
of an intermediate state, usually denoted by the 
letter C, prior to the formation of the precipitate 
zones. This compound is supposed to be produced 
continuously, with a certain rate proportional to 
the local concentration product of the reactants. 
It thus gives rise to a front of C particles that 
propagates in the system. Such processes can be 
tagged with the following symbolic reaction scheme 
: A + B — > C — > D. It is natural to suppose that 
this intermediate compound also diffuses inside the 
gel medium, with a diffusion coefficient D c . 

We recently considered two particular schemes of 
this type in a paper mainly devoted to a deeper 
understanding of the spacing law ]l0| . The first 
one, called the nucleation and growth theory, sup- 
poses that the intermediate compound undertakes 
a local nucleation process at the point where its 
concentration is greater or equal to a certain nucle- 
ation threshold c* . An extension of this theory to 
a case where counter-ions associated to the flow of 
the outer reactant play a catalytic role in the cre- 
ation of a new band was also examined under the 
denomination of induced sol-coagulation theory. 

As we noticed in Section II, the experiments we per- 
formed show unambiguously the presence of an interme- 
diate C compound. We will therefore base our theoreti- 
cal considerations on the class of models including the C 
species. 

B. Prediction of the intermediate-species class for 
the width law 

The determination of the width law for a particular 
model with intermediate species can be addressed at dif- 
ferent levels. 

A detailed description of the dynamics of band for- 
mation can be attempted. However, this route requires a 
precise knowledge of the microscopical mechanisms as ag- 
gregation, flocculation, etc. present in the system, which 
is clearly too ambitious a goal : very few informations 
are usually available about which of these mechanisms 
are present in a given experimental setup. 

Nevertheless, if one is only interested in general prop- 
erties of the bands system, as its width exponent a, a 
qualitative analysis based on general properties like con- 
servation laws should already provide a good approxima- 
tion for the width exponent. 

We are therefore in search of a general feature of all 
these intermediate species models that could allow us to 
compute a. Along this line, we found a general ar gum ent 
based on the shape and dynamics of the C front flQ|| . 

Let t n be the time at which band number n starts to 
appear. Over the interval [t n ,t n +i], the C profile adopts 
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a typical shape characterized by a diffusive tail at the 
right of its (growing) maximum, and a quasi-linearly de- 
creasing shape at the left, as shown in Figure 5. This fact 
can be easily verified by solving numerically the system 
of partial differential equations modeling the simplest 
intermediate-species theory producing Liesegang bands, 
namely : 



d t a(x,t) 


= D a dla( Xl t) 


— ka(x, t)b(x, t) 


(8) 


d t b(x, t) 


= D b d 2 x b(x,t) 


— ka(x, t)b(x, t) 


(9) 


d t c{x,t) 


= D c d^c(x, t) 


+ ka(x, t)b{x, t) — 


(10) 




Ki6(c-c*) 


c{x, t) — K2c(x, t)d{x, i) 




d t d(x, t) 


= Ki6(c-C*) 


c(x, t) + k,2c(x, t)d(x, t) 


(11) 



where k,K± and K2 are reaction rates (ki,K2 3> 1), and 
6 is the Heaviside step function. 

Note that in this simple basic model, we just require 
that C is rapidly going to zero wherever some precipitate 
already exists. This ingredient is aimed here at including 
crudely all the microscopic effects like aggregation. 

To discuss the width-law problem, we will balance the 
quantity of C diffusing out of the C front with the total 
amount of matter inside a band after its completion. 

Let's begin by summing up all the contributions to 
the quantity of C entering band n, once it is formed at 
.1 x n ■ They can be divided in two parts : 



c* 









\ C profile 



x n x n+1 position along 

the tube 



FIG. 5. Shape of the C profile just before band number n + 1 appears, x is the coordinate along the tube. 



Firstly, we can evaluate the quantity of C that flows 
into band number n during the time interval {t n ,t n+1 ). 
This quantity, denoted here Q n ,i, is equal to the time 
integral over the flux of C toward x n , evaluated from t n 
to t n +\- This flux is proportional to the slope of the C 
profile at the right of x n . Numerical estimations of this 
slope show that it scales with time as jl(| : 

d x C(x = x n + ,t n < t < t n+ i) 1 (12) 

v t — t n 

We therefore have 

D c d x C(x = X n + , t)dx — D c y/tn+l - t n 

(13) 

Invoking now the time and spacing laws (supposing that 
n is large enough), leads to 

Qn,l ~ D c y/x n+ i 2 - X n 2 — D c y/p 2 + 2px n S3 D c y^2px n 

(14) 

We go on then by examining the contribution Q ni 2 as- 
sociated to the flux into band n at times t > t n+ \. It 



turns out that once the next band appears at x = x n+ i, 
where the C profile has its maximum at t = t n +i, the 
left side of the front is trapped between two absorbing 
boundaries. At t = t n + , the total quantity of trapped C 
can be geometrically evaluated as (see figure 5): 

5C w -c*(x n+1 - x n ) (15) 

a finite fraction of which, say 77 (17 < 1), will diffuse into 
band n. Hence, using once again the spacing law, we can 
evaluate Q n ,2 as : 

Q n ,2 = ri5C w ~r]pc*x n (16) 
Finally, summing up the two contributions, we obtain : 
Qn = Qn,l + Qn,2 ~ x n (17) 

Note that typically, Q n ,2 is smaller than Q nt \ by four 
order of magnitude. 

The final step of our discussion deals with the relation 
between Q n and w n . 
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The basic model described by eq. (jq)-(|ll|) produces 
only infinitely thin bands. If one wants to generate bands 
with nonzero widths, it is necessary to supplement these 
equations with new terms accounting for growth mecha- 
nisms. Whatever the nature and details of these mech- 
anisms may be, one can make the following reasonable 
assumption : 

The terms one has to add to equations (|§|-|Tl|) in or- 
der to model growth lead, for t 3> i n +i, to a stationary 
profile of the precipitate concentration D inside the n th 
bands which adopts a natural scaling form: 



D = A n D(- 



Xt> "\ X < *~~, Xf 



(18) 



where A n is an amplitude which, in general, could depend 
on n. 

By equating the quantity Q n (the total amount of ma- 
terial that has nourished the band until its completion) 
with the total mass of precipitate inside the same band 
once formed, we have : 



dxD{- — — ) = jA n w n 



where gamma is a finite constant defined as 



1/2 



D(z)dz 



1/2 



To reproduce a relation of the form w n 
plitude A n should be of the form A n ~- 
hereby 

a + P = 1 



(19) 



(20) 



x", the am- 
!. We obtain 



(21) 



Thus the experimental value a = 1 corresponds (up to 
logarithmic corrections) to a constant amplitude A. This 
prediction could be tested experimentally by investigat- 
ing quantitatevely the internal structure of the bands. 

We can now return to the two theoretical predictions 
existing in the literature Jl],|llJ for the exponent a. In 
the first one, by Chopard et al., Lattice-Boltzmann cel- 
lular automata simulations including fluctuations were 
performed, in which microscopical mechanisms account- 
ing for aggregation were also included, in a phenomeno- 
logical way. They obtained bands patterns, which ex- 
hibited a width exponent varying between a — 0.5 and 
a = 0.6. Such a value, quite far from 1, can indeed be 
explained by the simple fact that their model allowed an 
unbounded quantity of precipitate to form at one lattice 
site. In this modeLthe amplitude A n increased with n 
and according to (|2l|), constrained the exponent a to be 
smaller than 1. 

In the second paper, Dee solved a system of mean-field 
equations similar to (||)-(11) in which the depletion term 
for the C front was build from physical considerations 
related to the nucleation and droplet growth theories. In 
this model, the distribution of precipitate inside a band 



as obtained from numerical computation seems to follow 
a scaling form with an amplitude approaching a constant 
after some transient regime. The width exponent that 
could be extracted from the set of (only) six bands pro- 
vided by the simulation is clearly compatible with the 
value a = 1. 



C. Summary and conclusion. 

We have reported pattern-forming experiments of 
Liesegang type, in which a careful measurement of the 
width of the bands leads to an exponent a close to 1. 
By considering the most basic theoretical model of the 
type related to the experiments we performed, we could 
explain, through a simple argumentation based on con- 
servation laws and macroscopic properties of the concen- 
tration profiles, how such a value close to unity can be 
understood, and which condition it implies on the shape 
of the precipitate profile inside bands. However, only a 
first principle derivation, based on the underlying micro- 
scopical mechanisms of coarsening could give a clear cut 
determination of a. 

We believe therefore that the microscopical theory of 
band formation and growth in Liesegang experiments 
constitutes an extremely rich and still unexplored sub- 
ject of investigation in itself, for which lots remains to be 
done. 
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